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Numerical calculations for the production of P-wave levels of Be quarkonium in 77 
collisions are performed in the leading O(asQfem) order of perturbation theory. The 
total cross-section of P-wave state production is about 10 % of that for the S-wave 
levels. The contribution of fragmentation component (6-1-6 diagrams) is low, and the 
basic contribution is determined by the recombination mechanism (8 Feynman diagrams). 
The gauge invariant term of the b —> Be fragmentation (6 diagrams) quite accurately 
reproduces the result of the fragmentation model, whereas there is a strong deviation of 
the c ^ Be fragmentation term from the predictions of the fragmentation model. 

I. INTRODUCTION 

At present, the production mechanism for the heavy quarkonium of mixed flavour. Be meson [1], is 
quite rehably predicted theoretically. The production of two pairs of heavy quarks bbcc is calculated in 
the framework of perturbative QCD, whereas the hadronization of the color-singlet be pair into the meson 
is described in a nonrelativistic potential model. 

In e"'"e" annihilation, the analysis of leading order approximation in the perturbation theory of QCD 
at M^/s <C 1 allows one to derive analytical expressions for differential cross-sections of the Be meson 
production being treated as the process of the b ^ Be fragmentation [2-5]. In the framework of the 
fragmentation mechanism, the S -wave level production dominates, and the P-wave state yield is only 
about 10 % with respect to the total number of produced mesons of the (be) family. However, as it 
was found in photon-photon and gluon-gluon collisions [6,7], there is an additional contribution of the 
recombination for the S'-wave states. Such a term changes spectra as well as the relative yields of 
pseudoscalar and vector states. The recombination mechanism can result in an enhancement of the 
P-wave level yield. 

In this paper we consider the exact calculation of complete set of the leading order diagrams in pertur- 
bation theory for the production of the (6c) system P-wave states in 77 collisions and make a comparison 
of the results with the fragmentation model approximation. 

II. CALCULATION TECHNIQUE 

The A^-^^' amplitude of the B^ meson production can be expressed through the amplitude of four free 
quarks production T^''' {pi, k{q)) and the orbital wave function of the Be meson, ^'■'^'^(q), in the meson 
rest frame as 

ASJ,. = f T^^^(p„fc(q)) . (vl.^'.(q))* .C//;^^, (1) 

where J and jz are the total spin of the meson and its projection on z axis in the Be rest frame, 
correspondingly; L and Iz are the orbital momentum and its projection; S and Sz are the sum of quark 
spins and its projection; C'^^[^ are the Clebsh-Gordon coefficients; pi are four-momenta of J5c, b and c, q 
is the three-momentum of b quark in the Be meson rest frame; /c(q) is the four-momentum, obtained from 
the four- momentum (0, q) by the Lorentz transformation from the Be rest frame to the system, where 
the calculation of T^'''{pi, fc(q)) is performed. Then, the four-momenta of 6 and c quarks, composing the 
Be meson, will be determined by the following formulae with the accuracy up to |q|^ terms 

Pb = TtPB. + ^(q): /o^ 
Pc=^PB.-fc(q), 

where rub and rUe are the quark masses, M = mt, + me, and Pb^ is the Be momentum. Let us note that 
for the P-wave states it is enough to take into account only terms, linear over q in eq.(0). 
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The product of spinors viUc, corresponding to the b and c quarks in the T'^"^ (pi, fc(q)) amphtude of 
eq.(|^), should be substituted by the projection operator 
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where F = 7^ for 5* = 0, or F = £*{Pb^, Sz) for S — 1, where £{Pb^, Sz) is the polarization vector for the 
spin-triplet state. 

For the sake of convenience, one can express the ViT) operator through the spinors of the following 
form 



u'ciPc - fc, ±) = ( 1 - 2^ U^ipc, ±), 



(4) 



where Vb{pb,i) and Uc(pc,i) are the spinors with the given projection of quark spin on z axis in the 
meson rest frame. Note, that the spinors in eq.(U) satisfy the Dirac equation for the antiquark with the 
momentum pb + k and mass rub or for the quark with the momentum Pc~ k and mass rric up to the linear 
order over k ( i.e. over q, too), correspondingly. 

One can easily show that the following equalities take place 
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^V{e*{P,+l)) + 0{k^). 
In the Be rest frame, the polarization vectors of the spin-triplet state have the form 

^'■^(-l) = 72 
e'^/(0) = (0,0,0,1), 
£'^^(+l) = -^(0,l,z,0). 



(6) 



In calculations the Dirac representation of 7-matrices is used and the following explicit form of the spinors 
is applied 
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For the P-wave states in eq.(|T]), the T'^"^ {pi, fc(q)) amplitude can be expanded into the Taylor series up 
to the terms linear over q. Then one gets 
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where R'p{Q) is the first derivative of the radial wave function at the origin, and has the following 
form 

'C°=afr' (9) 

~ \/2 V9«x ^dqy)^ 

where are the differential operators acting on T^^' {pi,k{q)) as the function of q = 

{<lx,qy,qz) at q = 0. 

As all considered matrix elements are calculated in the system distinct from the Be rest frame, the 
four-momentum fc(q) has been calculated by the following formulae 



k = q+( 1 _i^X:q,, (10) 



where v is the Be velocity in the system, where the calculations are performed. The matrix element 
T^^' (pi, k{q)) is computed, so that the four-momenta of b and c quarks are determined by eq.(|2|), taking 
into account eq.(p^). 

The first derivatives in eq.(P) are substituted by the following approximations 



9T^^'(R,fc(q)) , T^--(p,,fc(q^))-T^^-(p,,0) 
where A is some small value, and q-' have the following form 



(11) 



q" = (A,0,0), 

q^-(0,A,0), (12) 
= (0,0, A). 

With the chosen values of quark masses and interaction energies, the increment value A = 10^^ GeV 
has provided the stability of 4-5 meaning digits in the squared matrix elements summed over for all 
P-wave states with the given value of J and S, when one has performed the Lorentz transformations 
along the beam axis or the rotation around the same axis. 

One has to note that because of such transformations, the new vectors fc(q-') do not correspond to 
the transformed old vectors. Therefore, the applied test is not only a check of the correct typing of the 
T^^' {pi, k{q)) amplitude, but it is also the check of correct choice of the phases in eq.(^). 

The matrix element A^"^^' squared, which is calculated by the method described above, must be 
summed over jz as well as the spin states of free b and c quarks. It also must be averaged over spin 
projections of initial particles. 

The phase space integration has been made by the Monte Carlo method of RAMBO program [8] . 

III. DISCUSSION OF RESULTS 

To check the numerical way of the amplitude calculation for the P-wave level production of Pc, we 
have considered the production of these states in e+e" annihilation, where the analytical expressions 
for the differential cross-sections were derived in the M^/s <C 1 limit [5]. Those expressions define the 
functions of the b ^ Bc{L = 1) fragmentation. It should be noted that in the accurate consideration of 
the fragmentation mechanism for the Pc(P = 1) production in e+e~ ^ ^* ^ Be + X , one can see that 
in addition to the b fragmentation one has to account for the c quark fragmentation into Be- Moreover, 
the most significant role is played by the c fragmentation into '^Pq state. Therefore, to compare with the 
numerical results we use the analytical expressions accounting for both b and c fragmentation. As one 
can see in Fig.l, the distributions of ^Pi, "^Pq, '^Pi, '^P2 level R production, da/dz (z = 2|PbJ/y^ with 



^In the fee system the quark spin-dependent corrections to the heavy quarkonium potential lead to the mixing 
of ^Pi and ^Pi levels [9]. In the fragmentation model this mixing results in redefinition of corresponding frag- 
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being the three-momentum of Be meson in the c.m.s.), calculated numerically and given analytically, 
coincide with each other at the same set of parameters, which are give below 



aem = 1/128, 
as = 0.2, 

rrib = 5.0 GcV, (13) 
nic = 1.7 GeV, 
|i?^(0)|2 = 0.201 GeV^ 

Thus, the performed verification convinces us that the calculation method used is quite accurate. 

Let us consider the P-wave level production of Be in 77 collisions. The cross-section of the P-wave 
level production at various energies of interacting photons are presented in Tab.l and Fig. 2. One can see 
in Fig. 2 that in the region of interest the energy dependence of the summed cross-section at the chosen 
TOfc and rric values is quite accurately described by the following approximation, shown as solid line in 
Fig.2, 

/ \2\^'''/ \1 32 

2{ini, + mc)\ \ f 2{mi, + rric) 



...(..r) = 130 . 1^1 - ) j • ) (14) 

The summed fragmentation contribution, obtained as the product of the 77 bb and 77 cc cross- 
sections and the corresponding probabilities of the fragmentation (5.34- 10~^ for the b quark fragmentation 
and 1.58-10"^ for the c quark one, respectively), is shown in the same figure. One can see in Fig.2 that the 
fragmentation contribution, evaluated in the model, is small at high energies, where the application of this 
model could be sound for the b fragmentation, at least. The fragmentation model overestimates the exact 
result obtained over the complete set of leading order diagrams at low energies close to i/i < 30 GeV. 
This overcstimation has a simple explanation and it is related with incorrect evaluation of phase space, 
since in the fragmentation model, one uses the two-particle phase space instead of the three-particle one 
in the exact calculations. For the correct study of the fragmentation mechanism the M^/s ^ 1 condition 
is necessary to be satisfied. Therefore, in what follows, we will restrict ourselves by the consideration of 
differential distributions at = 100 GeV, where the mentioned condition is certainly valid. 

The da/dz distribution for ^Pi state production cross-section in 77 collisions is shown as the solid 
line histogram in Fig. 3a. The dashed line histogram in the same figure denotes the gauge invariant 
contribution of six diagrams, where the cc pair is emitted from the 6 or 6 quark line, i.e. the b quark 
fragmentation into Be meson takes place. This exactly calculated term is compared with the prediction 
of fragmentation model, considered in [5] and presented as the smooth dashed curve. One can see that 
the exact result for b ^ Be is quite accurately described by analytical expression of fragmentation 
model, whereas the c fragmentation diagrams contribution, shown as the dotted line histogram, is larger 
than the predictions of fragmentation model, so that in absolute values, it is larger than the b quark 
fragmentation term. Remember, that the analogous picture takes place also for the photonic production 
of S'-wave states of Be meson [6]. One can sec in Fig. 3a, that as well as in the production of 5- wave 
levels, the fragmentation contribution does not dominate, and the main contribution is determined by 
the recombination diagrams. However, for the correct study of the fragmentation mechanism in photon- 
photon collisions, one must consider the spectra not over the total energy of the meson, but over its 
transverse momentum, giving an additional scale of energy, so that one can expect the factorization of 
the heavy quark hard production at large Pt ^ Mb^ and the forthcoming fragmentation, where the 
particle virtualities are of the order of the quark masses. 

The cross-section distributions over the transverse momentum of ^Pi level of the Be meson are shown in 
Fig. 3b. The solid line histogram denotes the result of calculations over the complete set of diagrams, the 
dashed histogram and curve show the contribution by the b fragmentation diagrams and the prediction 
of fragmentation model, correspondingly, and the c fragmentation diagram contribution (the dotted 



mentation functions, i.e. in the introduction of additional functions for the 1^ and l"*"' states [5]. However, first, 
the mixing effect does not influence the consideration of the physical mechanism of photonic production of the 
P-wave states, and second, in the studied approach the '^Pj-state masses are degenerated over J, so that the 
isolation of the S = and S = 1 components in the J = 1 state is rather conventional. That is why we restrict 
ourselves by the consideration of distributions for the ^ Pj states. 
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histogram) is compared with the corresponding prediction of fragmentation model (the dotted curve) at 
y/s = 100 GeV. 

As well as in the S-wave state production, the fragmentation mechanism does not dominate even at 
large Pt- 

The analogous distributions for the '^Pq, ^Pi, ^P2 states are presented in Figs.4a-6b. In Figs. 5-6, where 
the ■^Pi, ^P2 spectra are shown, one can draw the conclusions, which repeat the statements concerning 
the ^Pi production. The picture of the '^Pq meson production slightly differs from the general case. As 
one can see in Fig. 4b, the b quark fragmentation at large transverse momenta plays a much greater role 
than that in the production of other P-wave levels, while, contrary, the diagrams corresponding to the c 
fragmentation, are less essential than in other cases. 



IV. CONCLUSION 



In this paper we have performed the numerical calculation for the P-wave level production of Be mesons 
in 77 collisions in the leading 0(alal^) order of the perturbation theory. From the theoretical point 
of view, the consideration of 77 collisions in this respect is of special interest, since, on the one hand, 
in the photonic production of Be as well as in e^e" annihilation, one can isolate the gauge invariant 
set of fragmentational diagrams, and, on the other hand, there is also the contribution of recombination 
type diagrams, which are essential in the consideration of a more complicated case of the hadronic Be 
production [7]. 

The performed calculations show that: 

1. The total cross-section of P-wave state production of Be is about 10 % in respect to the production 
of 5-wave levels, as it takes place in e+e~ annihilation. 

2. As well as in the S'-wave level production, the recombination mechanism dominates, and the frag- 
mentation one is small. 

3. The h fragmentation diagram contribution is quite accurately described by the fragmentation func- 
tion at large Pr as well as low one, whereas for the c fragmentation diagrams, the fragmentational 
picture is broken completely. 
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The work of A.V. Berezhnoy has been made possible by a fellowship of INTAS Grant 93-2492 and one 
of International Soros Science Education Program Grant A1377 and is carried out within the research 
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TABLE I. The dependence of production cross-sections for different P-wave states of Be on the total energy 
of colliding photons. (The calculation errors are shown in parenthesis). 



Vi, GoV 


o-ipj, fb 


fT3po, fb 


(73 , fb 


173 pj, fb 


15. 


1.064(2) 


0.1329(3) 


0.08707(18) 


0.894(2) 


20. 


5.703(11) 


0.968(2) 


1.147(2) 


8.098(15) 


30. 


7.56(3) 


1.409(5) 


2.84(1) 


12.11(5) 


40. 


6.72(4) 


1.269(9) 


3.063(18) 


11.01(7) 


60. 


4.68(5) 


0.882(12) 


2.42(3) 


7.7(1) 


80. 


3.32(5) 


0.624(15) 


1.80(3) 


5.48(12) 


100. 


2.58(3) 


0.493(9) 


1.42(2) 


4.32(7) 
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FIGURE CAPTIONS 



Fig, 



1. The Be production cross-section distributions over z in the e+e 



7* 



Be + X process are 



presented as the histograms in comparison with the fragmentation model predictions shown as the 
smooth curves, at y/s = 100 GeV for the following Be states: ^Pi (solid line), '^Pq (dashed line), 
^Pi (dotted line), (dotted-dashed line). 

Fig. 2. The summed cross-section dependence of P-wave states on the energy of interacting photons is 
marked by (•) in comparison with the prediction of b and c q uark fragmentation model (dashed 
curve). The solid line corresponds to the approximation of (|l4|). 



Fig. 3. a. The z distributions, corresponding to the photonic production of ^Pi state at the interaction 



energy 100 GeV. The total result is presented by the solid line histogram, the b fragmentation 
diagrams contribution (dashed histogram) is compared with the prediction of fragmentation 
model (dashed curve), the c fragmentation diagrams term is denoted as the dotted histogram 
in comparison with the fragmentation model result (dotted curve). 

b. The transverse momentum distributions, corresponding to the photonic production of ^Pi 
state at the interaction energy 100 GeV. The notations of different contributions are the same 
as in Fig. 3a. 




Fig. 4a, b. The ■^Pq spectra, denoted as in Fig. 3a, b. 
Fig. 5a, b. The '^Pi spectra, denoted as in Fig. 3a, b. 
Fig. 6a, b. The '^P2 spectra, denoted as in Fig. 3a, b. 
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Fig.l 
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